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ABSTRACT 

An exact solution is presented to a model that mimics the crowding effect in financial 
markets which arises when groups of agents share information. We show that the size 
distribution of groups of agents has a power law tail with an exponential cut-off. As the 
size of these groups determines the supply and demand balance, this implies heavy tails 
in the distribution of price variation. The moments of the distribution are calculated, 
as well as the kurtosis. We find that the kurtosis is large for all model parameter values 
and that the model is not self-organizing. 
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1. Introduction 

Empirical analyses of financial price-data show that short time variations of the 
price of different assets deviate from the Gaussian distribution!!], which would be 
expected if agents were trading independently. Anomalously large oscillations are 
far more probable than predicted by a Gaussian distribution and are distributed 
according to an exponentially truncated power lawo. 

Physically motivated models used to describe this phenomena assume that it is 
caused by cooperative behaviour of the tradersa'cl. These models make use of the 
properties of system in a critical state, where the quantity describing the range of the 
interactions between agents is diverging)^. Generally, in a critical state, the elements 
that make up the system behave synchronously and a microscopic formulation is no 
longer required to describe kinetics of the whole systenJj. 

Having made this observation, Cont and BouchaudB proposed a model of ran- 
domly connected agents to describe crowding in financial markets. Connected agents 
are agents sharing the same information and making the same decision between buy- 
ing, selling or doing nothing. In their original model, any pair of agents is connected 



with a probability p, corresponding to bond percolation in an infinite dimensional 
spaceo or cquivalently, to agents being vertices of a random graphO. To drive the 
system into a critical state, the probability p that two agents are connected has to 
be inversely proportional to the number of agents, Nq. Assuming that p is pro- 
portional to 1/Nq, they find a power law distribution for the cluster size, with an 
exponential cut-off, in reasonable agreement with the empirical data. When p is 
exactly equal to I/Nq, the exponential correction disappears. 

We present in this paper an exact solution to an extension of the Cont-Bouchaud 
model!, introduced by Egufluz and Zimmermanni. Instead of supposing that two 
agents are connected with a probability p, they considered that the network of 
agents is building itself and changing dynamically as decisions are made by the 
agents. Egufluz and Zimmermannel investigated numerically the simplest case where 
the network of agents is built randomly. As expected, they obtained numerically 
the same distribution for the returns as Cont and Bouchaudtl. 

In Sec. 2, we present the model of Egufluz and Zimmermanni and find an exact 
analytical expression for the size distribution of the clusters of agents. A cluster 
of agents is a set of agents making the same decision, between buying, selling and 
doing nothing. We a^ume that the traded amount is proportional to the number 
of agents in a clustera. Then, the distribution of returns, which is the distribution 
of the relative difference between the number of buyers and the number of sellers, is 
a simple function of the cluster size distribution. In Sec. 3, a generalisation of the 
model is proposed. Although the cluster size distribution is not obtained exactly in 
this case, its asympotic behaviour is determined and a method to obtain its various 
moments is explained. 

2. The Model 

Consider a set of No agents in which connected agents are able to exchange 
information. An agent can be connected to any of the No — 1 other agents. When 
it is possible to go from one agent to another following a path of connected agents, 
these two agents are said to belong to the same cluster, that is, a cluster is a set of 
agents that can exchange information, either directly or through the intermediary 
of other agents. All the agents are isolated at the beginning of the simulation. At 
each time step, an agent is selected at random. With a probability a, this agent 
decides to buy or sell. In this case, all the agents in the same cluster follow the lead 
agent in buying or selling and the cluster is broken up into isolated agents. With a 
probability 1 — a, he decides to do nothing. In this case, two agents are selected at 
random and connected to each other. 

When an agent decides to make a transaction, buying or selling, he triggers the 
action of all the agents belonging to the same cluster. If the selected agent is buying 
or selling, all the agents belonging to the same cluster are buying or selling at the 
same time, respectively. Hence, the variation of the supply and demand balance, 
that is, the return, is a function of the size distribution of the clusters of agents. 



In this model, the return R at a given time is defined as the ratio of the agents 
buying or selling at this time. If the agents are buying, the return is positive and 
it is negative if the agents are selling. Hence, if n s is the number of clusters of size 
s, the probability to have a return of size s/Nq is given by sn s /No. 
The number n s of clusters of size s > 1 evolves like 



where N is the total number of agents. Note that one time step in this formulation 
corresponds to one attempted update per agent in the numerical simulation. The 
first term on the right hand side of Eq. (|l|) describes the fragmentation of a cluster 
of size s in case of a transaction involving s agents. The second term describes the 
creation of a new cluster of size s by coagulation of two clusters of size r and s — r, 
with r < s. The last term on the right hand side describes the disappearance of a 
cluster of size s by coagulation with another cluster. The number of clusters of size 
1, rti obeys 



The first term on the right hand side describes the appearance of s clusters of 
size 1 by fragmentation of a cluster of size s while the second term describes the 
disappearance of clusters of size 1 by coagulation with another cluster. 

When the system reaches a stationary state, the size distribution of the clusters 
obeys 
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for s > 1 and 
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Introducing the generating function 



oo 




(•5) 



r=2 



we find 




giving a size distribution for the clusters of 



Using Stirling's formula to expand this for large s, we have 

Hence, the model presents a power law distribution of exponent a = 5/2 for the size 
of the clusters, with an exponential cut-off, as in the model of Cont and Bouchaudtl. 
For a = 0, the exponential cut-off vanishes. 

The properties of the cluster size distribution n s can be calculated exactly. The 
number of clusters per site, Mo, is given by 

M = 2-fc^ln(2-a). (9) 
(1-a) 

This is easily obtained by integrating g(u>). It is interesting to note that, even if 
the probability of a transaction is 0, the system does not become a single cluster 
containing all the agents. In fact, a — > is the only value of a that gives a power 
law for the cluster size distribution. The other moments of the distribution, Mi, 
defined as 

^ oo 

can be obtained by differentiation of the generating function. We find that M\ = l 
and M2 = 1 1 a. As the probability of a return of size s/No is given by sn s /No, the 
variance of the distribution of returns, er 2 , is 

a 2 = A/ 3 -M|= (1 ' fl) l 2 ' a) . (II) 
Finally, the kurtosis of the distribution of returns is given by 



a 4 - 18a 3 + 78a 2 - 120a + 60 
o(l - a)(2 - a) 

where /i4 is the fourth central moment of the distribution of returns. The kurtosis 
k is shown as a function of a in Fig. 1. It diverges at a = and a = 1, while 
being asymmetric around a = 1/2. k is not small for any value of a, achieving a 
minimum value of k ~ 26.07 for a ~ 0.846. This shows that, even if the cluster size 
distribution is not a power law for a > 0, it has heavy tails for every value of a. 
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Fig. 1. Analytical result for the kurtosis n of the return distribution as a function of a for m = 2. 
Note that the ordinate is a logarithmic scale. 



3. Generalisation 

The generalisation wc introduce is different from the model of the previous 
section in that when an agent decides with a probability 1 — a to do nothing, m 
agents are selected at random and connected together. Then, the model of the 
previous section corresponds to m = 2. 

If the simulation starts with isolated agents, the possible sizes of the clusters 
are given by (m — l)i + 1, where i is a non-negative integer. It is easy to show that 
the combination of m such clusters also gives such a cluster. Hence, these are the 
only possible sizes. Excluding all the other cluster sizes, the time evolution of the 
number of clusters of size s > 1 is given by 
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with a particular equation for s = 1 



m — 1 

(15) 

When the system has reached a stationary state, the cluster size distribution is the 
solution of 
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where the Kronecker delta S a .f3 = 1 when a = (3 and when a p. 

Considering m = 2, the equation giving the stationary distribution can also be 
written 
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by inserting the expression for rn r into the expression of sn s . Consequently, when 
s is large, the stationary distributions for to = 2 and to = 3 display the same 
dependence on s. Introducing iteratively the expression for rn r in sn s , it is easy to 
show that the stationary distribution displays the same dependence on s for every 
value of to when s is large. From the previous section, the stationary distribution 
is a power law with an exponential cut-off. The exponent is a — 5/2, independent 
of to. Considering sn s for the first values of s, as obtained by Eq. (|l6|), the cluster 
size distribution also displays a dependence like 

/ 1 \ \ (s — l)/(m — 1) / \ r 

(1-a) y I — I 
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By iteration, sn s follows this relation for all values of s, giving an exponential 
correction to the power law, as expected. Fig. 2 compares numerical results for sn s 
for to = 2 and to = 3, using Eq. (|l6| ) iteratively, with a = 0. An exponent of 5/2 is 
obtained in both cases. Note that in direct simulations of the model, it is difficult 
to appreciate the departure from the power law from the numerical simulations of 
the model for a < 0.5, as the largest values of s are restricted. A power law fit to 
the results from direct numerical simulations of the model for to = 2, No = 10 4 and 
a = 0.01 gives an exponent!] a ~ 2.7. Hence, it is necessary to consider much larger 
systems to make a correct estimate of a. 

Introducing the same generating function g(u>), Eq. we obtain the self- 
consistent equation 
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Fig. 2. Numerical results for the cluster size distribution n a for m = 2 (upper dots) and m = 3 
(lower dots) obtained by iterating Eq. (|l(]) with a = 0. For in = 2, n s has been rescaled by a 
factor of 100 to separate the two sets of results. 
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a + m{l — a) 

By differentiating Eq. (^0|), all the moments of g(co) can be calculated exactly. For 
instance, 
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4. Conclusions 

We have presented the exact solution to a model that mimics the herding effect in 
financial market as well as the information transmission. A power law distribution 
with an exponential cut-off is found for the size distribution of the clusters of agents 
and hence for the distribution of returns. The model is not self-organized as the 
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parameter a describing the probability of a transaction has to be tuned to a low 
value to ensure the appearance of transactions of all sizes. Except for the small 
values of a, the system does not display a power law distribution for the large 
transactions. Nevertheless, we computed the kurtosis to show that the model has 
heavy tails for any value of a. Hence, the model deviates from Gaussian behaviour 
and the assumption of independent agents. 

A simple generalisation of the model was also investigated analytically, consid- 
ering the number of agents connected at a time step as a free parameter to. The 
large s dependence of the distribution of the clusters of agents is determined. The 
cluster size distribution is found to be independent of m in the large s limit. A 
method to calculate exactly any moment of the distribution was presented. 

To incorporate more realistic features, several extensions of this model are pos- 
sible, such as a variable to or a different process of cluster fragmentation, allowing 
some agents to remain connected after a transaction. However, the modifications to 
the equations seems to be irrelevant for large values of s. Consequently, the same 
asymptotic behaviour with a power law of exponent a = 5/2 for the cluster size 
distribution is expected. 
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